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Abstract 

It is proved that the measurement of the acoustic pressure on the ear membrane 
allows one to determine the shape of the ear canal uniquely. 

1 Introduction 

Consider a bounded domain D C M. n , n = 3, with a Lipschitz boundary S. Let F be an 
open subset on S, a membrane, G = S \ F, T = dF, and N is the outer unit normal 
to S. The domain D models the ear canal. The acoustic pressure u is given on the 
membrane F, u = f ^ on F. This pressure is zero on G, and the normal component 
of the velocity, '■= h, is measured on F. The problem we are concerned with is: given 
the datum {f,h}, with f ^ 0, for a single f, can one recover uniquely the shape of D, 
the ear canal? 

Let us introduce now the corresponding mathematical formulation of the above inverse 
problem. Consider the problem: 

V 2 m + k 2 u = in D, u = f on F, u — 0onG. (1.1) 

We assume that k 2 is not a Dirichlet eigenvalue of the Laplacian in D. This assumption 
will be removed later. If this assumption holds, then the solution to problem (1.1) is 
unique. Thus, its normal derivative, un '■= h on F, is uniquely determined. Suppose one 
can measure h on F for a single / G C 1 (F), f ^ 0. 

The inverse problem (IP) we are interested in can now be formulated: 

Does this datum determine G uniquely? 

Thus, we assume that F, f and h are known, that k 2 is not a Dirichlet eigenvalue of 
the Laplacian in D, and we want to determine the unknown part G of the boundary S. 
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Let A be the smallest eigenvalue of the Dirichlet Laplacian L in D. Let us assume 
that 

A > k 2 . (1.2) 

Then problem (1.1) is uniquely solvable. Assumption (1.2) in our problem is practically 
not a serious restriction, because the wavelength in our experiment can be chosen as we 
wish. Since the upper bound on the width d of the ear canal is known, and since 

A > 1. (1.3) 

one can choose k 2 < -4 to satisfy assumption (1.2). A proof of the estimate (1.3) is given 
at the end of this note. 

We discuss the Dirichlet condition but a similar argument is applicable to the Neu- 
mann and Robin boundary conditions. Boundary- value problems and scattering problems 
in rough domains were studied in [1]. 

Our basic result is the following theorem: 

Theorem 1. If (1.2) holds then the above data determine G uniquely. 

Remark 1. If k 2 is an eigenvalue of the Dirichlet Laplacian L in D, and m(k) is 
the total multiplicity of the spectrum of L on the semiaxis A < k 2 , then G is uniquely 
defined by the data {fj, /ij}i<j< m (fc)+i, where {fj}i<j< m {k)+i is an arbitrary fixed linearly 
independent system of functions in C(F). 

In Section 2 proofs are given, in Section 3 a numerical approach to computing G is 
discussed very briefly, and in Section 4 conclusions are formulated. 



2 Proofs. 

Proof of Theorem 1. 

Suppose that there are two surfaces G\ and G2, which generate the same data, that 
is, given / on F they generate the same function h on F. Let Di,u\ and D 2 ,u 2 be the 
corresponding domains and solutions to (1.1). Denote w := U\ — u 2 , D 12 := D\ D Di, 
D 12 := D 1 U D 2 , D 3 := D 1 \ D 12 , £> 4 := D 2 \ D 12 . Note that w = w N = on F, since 
the data / and h are the same by our assumption. Consequently, one has: 

V 2 w + k 2 w = in D 12 , w = w N = on F. (2.1) 

By the uniqueness of the solution to the Cauchy problem for elliptic equations, one 
concludes that w = in D 12 . Thus, ui = u 2 = on dD 12 , and u\ = on dD 3 . Therefore 

VV + k 2 u x = in D 3 , Ul = on dD 3 . (2.2) 

Since D 3 C D, it follows that A(D 3 ) > A(D) > k 2 . Therefore k 2 is not a Dirichlet 
eigenvalue of the Laplacian in D 3 , so u x = in D 3 , and, by the unique continuation 
property for solutions of the homogeneous Helmholtz equations one concludes that u± — 
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in D]_. In particular, U\ = on F, which is a contradiction, since U\ = f ^ on F by 
the assumption. Theorem 1 is proved. □ 

Proof of Remark 1. Suppose that k 2 > is arbitrarily fixed, and the data are 
{/?'; hj}i<j<m(k)- Using the same argument as in the proof of Theorem 1, one arrives 
at the conclusion (2.2) with Uj t \ in place of ui, where solves (1.1) with / = fj, 
1 < j < m(k) + 1. Since the total multiplicity of the spectrum of the Dirichlet Laplacian 
in D is not more that m(k), one concludes that D\ = D 2 . Remark 1 is proved. □ 

Proof of estimate (1-3). 

For convenience of the reader and to make the presentation self-contained, we give 
a short proof of inequality (1.3), although this inequality is known. Let a be a unit 
vector, and d(a) be the width of D in the direction a, that is, the distance between two 
planes, tangent to the boundary S of D and perpendicular to the vector a, so that D 
lies between these two planes. Let 

d : = mind(a) > 0. 

a 

By the variational definition of A one has: A = min j D | Vu\ 2 dx, where the minimization is 
taken over all u G H 1 , vanishing on S and normalized, | \u\ \l 2 (d) — 1- Denote s := x±, y : = 
(x2, xs), and choose the direction of x\— axis along the direction a, which minimizes d(a), 
so that the width of D in the direction of axis X\ equals d. Then one has: 

u(s,y) = / u t (t,y)dt, 

J a 

,y)?< \u t (t,y)\ 2 dt(s-a)<d \u t (t,y)\ 2 dt, 

J a J a 

where s = a and s = b are the equations of the two tangent to S planes, the distance 
between them is d = b — a, and D is located between these planes. 

Denote by F s the crossection of D by the plane x\ = s, a < s < b. Integrating the 
last inequality with respect to y over F s , and then with respect to s between a and b, 
one gets: 



so 

\u(s 



\u\ 



122(D) <d ii v«ir, 



which implies inequality (1.3). □ 



3 A discussion of the numerical aspects of the prob- 
lem 

Let us discuss very briefly how to calculate G numerically, given the data {f,h}. We 
have to calculate three scalar functions which give a parametric equation of G, of the 
form Xj = <fij(t, v), j = 1, 2, 3, t, v e (0, 1) x (0, 1) are parameters, the three functions 4>j 
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are unknown. They are to be found from the three conditions: u — f on F, = h on 
F, and u = on G. One writes a representation for u in D by Green's formula 

u= [g(x,s)h(s) - g N (x,s)f(s)]ds+ / g(x, s)u N (s)ds := U(x) + V(x), 
Jf Jg 

where g = f^zjr, the function «jv on G is not known, ds is an element of the surface 
area, the functions <pj, j = 1,2,3, are not known, the function U(x) is known, and 
V(x) is a single-layer potential with the unknown density un '■= H on G. To determine 
the four unknown functions <f>j,j = 1,2,3, and H, one has three boundary equations: 
u = f on F, Ujsi = h on F, u = on G, and the fourth equation for H on G one 
obtains from the formula H = Un + \AH + \H for the normal derivative of the single- 
layer potential on G, which yields an integral equation for H: H = AH + Un- Here 
AH := 2 f G ^jv(S) s')H(s')ds'. One expresses (is in terms of the functions <pj, and get four 
nonlinear integral equations for four unknown functions <pj, j = 1,2,3, and H. These 
equations do have a solution if the data {/, h} are exact. The solution can be obtained 
by a Newton-type method. The problem is ill-posed, because small variations of h lead 
to a function which may be not a normal derivative on F of the solution to problem (1.1) 
corresponding to the given function /. Thus, one has to use a regularization for solving 
the above integral equations numerically. 

4 Conclusion 

The basic result of this note is the proof of the following statement: if one applies 
some pressure to the eardrum and measures the corresponding normal component of the 
velocity of this drum, then one can uniquely determine from these data the shape of the 
ear canal. This note is of theoretical nature, but hopefully it may lead to a progress in 
the construction of better hearing devices. 
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